Regular and chaotic transport of discrete solitons in asymmetric potentials 
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Ratchet dynamics of topological solitons of the forced and damped discrete double sine- Gordon 
system are studied. Directed transport occurring both in regular and in chaotic regions of the 
phase space and its dependence on damping, amplitude and frequency of the driving, asymmetry 
parameter, coupling constant, has been extensively investigated. We show that the passage from 
ratchet phase-locked regime to chaotic ratchets occurs via a period doubling route to chaos and 
that, quite surprisingly, pinned states can exist inside phase-locking and chaotic transport regions 
for intermediate values of the coupling constant. The possibility to control chaotic discrete soliton 
ratchets by means of both small subharmonic signals and more general periodic drivings, has also 
been investigated. 
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I. INTRODUCTION 

Ratchet dynamics in nonlinear systems are presently 
attracting a great deal of interest due to their relevance in 
several fields, including molecular motors and Josephson 
junctions [H[2]. For point particles ratchets the directed 
motion achieved in presence of damping and determinis- 
tic forces of zero average, is a direct consequence of the 
breaking of the spatio-temporal symmetries in the system 
[3HS]- The resulting net motion of the particle is usually 
phase-locked to the external driving [5] and the deter- 
ministic ratchet dynamics is quite robust to exist also in 
presence of noise. 

Ratchet motion occurring in infinite-dimensional con- 
tinuous nonlinear systems with soliton-like solutions 
(soliton ratchets) have been extensively investigated in 
the past years. This has been done both in the case of 
spatially asymmetric potentials and single harmonic driv- 
ing [21 [TO] and for symmetric potentials and temporarily 
asymmetric bihharmonic forces [51 111H13) . Experimen- 
tal implementations of soliton ratchets in long Josephson 
junctions were also done in |14j using time-asymmetric 
biharmonic currents, and in [15] and |16) using asym- 
metric magnetic helds or spatially asymmetric currents. 

The mechanism underlying continuous soliton ratch- 
ets, also known as internal mode mechanism, is ascribed 
to the coupling between the soliton internal degree of 
freedom {internal mode) and the translational mode of 
the soliton [TU]. This coupling is induced mainly in the 
presence of damping and allows an asymmetric transfer 
of energy from the external driving to the internal mode 
and to the soliton center of mass which generates the di- 
rected motion. Thus, opposite to point particle ratchets 
for which the action of the force is always direct on the 
center of mass, in soliton ratchets the force acts indirectly 
on the center of mass and involves the vibrational mode 
of the solitons in an essential manner (a direct action 
may exist in the overdamped limit where the phenom- 
ena becomes negligibly small [9]). This is quite natural if 
one think that the soliton is not a rigid object but it can 



vibrate while moving, just like a droplet of liquid does 
when dropping. This internal mode mechanism has been 
confirmed to be of general validity for continuous dissi- 
pative soliton systems [SJ [TUJ [TTH20] . The possibility of 
soliton ratchets in continuous Hamiltonian systems has 
been recently demonstrated for the case of the nonlinear 
Schrodinger equation with asymmetric potentials rapidly 
flashing in time [21] . 

In contrast to the continuous case, discrete soliton 
ratchets (DSR) in nonlinear lattices have been much 
less investigated. Detailed studies exist only for ratchet- 
like behaviors in a dissipative discrete sine-Gordon chain 
driven by a biharmonic force which breaks the time sym- 
metry [22] [23] , and for the case of a conservative discrete 
nonlinear Schrodinger equation describing a ID array of 
coupled waveguide optical resonators [24] . In both cases 
the interplay between discreteness and nonlinearity in- 
troduce new features in the phenomenon such as: ex- 
istence of depinning thresholds for transport, stair-case 
like dependence of the mean drift velocity on system pa- 
rameters, diffusive and even intermittent ratchet-like dy- 
namics [22] . The possibility to enhance DSR transport 
by introducing phase disorder into the asymmetric bihar- 
monic periodic driving was also considered [25] . 

It is interesting, however, to investigate the phe- 
nomenon of DSR in the case in which the potential is 
asymmetric in the field variable and the system is driven 
by harmonic forces of zero mean. In this case the dis- 
crete solitons have an internal degree of freedom even in 
absence of damping and driving, a fact which could lead 
to qualitative different behaviors. 

The aim of this paper is twofold. From one side, we ex- 
tend previous investigations on DSR to the case of asym- 
metric potentials. For this we consider the case of a dou- 
ble sine-Gordon chain with potential function given by 

V{x) = C - cos(27rx) + ^ sin(47rx), A G [0, 1], 

(1) 

where A is an asymmetry parameter and C is a con- 
stant that fixes the zero of the potential. Therefore the 
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system we are going to consider belongs to the Frenkel- 
Kontorova family |25j and can be visualized as a chain of 
mechanical double pendula connected by a gear of ratio 
1/2 with nearest-neighbor interactions and with the peri- 
odic on-site double-sine potential in ([!]) due to the gravity 
[26j . The dependence of the transport on the damping 





FIG. 1: (Color online), (a) Comparison between the symmet- 
ric potential for A = (continuous line) and an asymmetric 
double sine-Gordon potential (dashed line) for A = 0.5. (b) 
Stationary discrete kink solutions of the unperturbed double 
sine-Gordon chain corresponding to A = (circles joined by a 
continuous line) and A = 0.5 (squares joined by a continuous 
line). The discrete spatial derivatives of the kink profiles are 
also shown in (c) with corresponding symbols. The value of 
the coupling constant k is fixed to 1. 



constant, amplitude and frequency of the driving force, 
coupling constant and asymmetry parameter, is investi- 
gated in detail. As a result we show that in contrast 
to DSR driven by asymmetric forces of zero mean [2"2"] . 
kink transport in the present system can be quite effec- 
tive not only in correspondence of regular (phase-locked) 
orbits,but also on chaotic trajectories. 

/From the other side, we shall investigate the DSR 
chaotic transport in some detail both with respect to 
the onset of chaos and with respect to the possibility of 
controlling the chaotic ratchet trajectories using subhar- 
monic signals which do not break the time symmetry in 



the system. To this regard, we introduce a suitable map 
in the space of the kink coordinates by means of which 
we construct bifurcation patterns. Using this approach, 
we show that the chaos onset occurs via period doubling 
and that the chaotic region can be strongly reduced in 
favor of periodic phase-locked orbits in presence of small 
subharmonic signals. The complete suppression of chaos 
is also investigated by means of elliptic driving functions 
with modulus taken as a free parameter to achieve op- 
timal control. For the regions of parameters we have 
investigated, we show that the full chaos suppression can 
be implemented in any practical context by means of two 
properly chosen harmonic signals. 

The paper is organized as follows. In Section II we in- 
troduce the model equations and discuss the main prop- 
erties of the static kink solutions, including their Peirls- 
Nabarro barrier and the small oscillation problems. In 
section III we investigate the dependence of the discrete 
ratchet dynamics, both phase-locked and chaotic, on the 
system parameters. In Section IV we concentrate on the 
ratchet dynamics associated to chaotic orbits and on their 
stabilization by means of small subharmonic modulations 
of the harmonic driving. We study the onset of chaos by 
means of a suitable bifurcation map and show that chaos 
occurs via period doubling. Finally, in Section V we sum- 
marize the main conclusions of the paper. 



II. MODEL EQUATION 

The chain of double pendula associated to potential 
([l]) is described in a dimensionless form by the equation 

i„+ai n +U'(x„)+K(2a;„-a; n+ i-x„_i) = J+F(t)+£ n (t), 

(2) 

with n € [-N, N], x n being the displacement of the nth 
pendulum with respect to its equilibrium position, 7 a 
constant bias term, a the damping coefficient, k the cou- 
pling constant which measures the discreteness of the lat- 
tice, 



F(t) 



2tt 



cos(ujt + 9) 



(3) 



an external symmetric in time periodic driving force of 
strength e, frequency uj and £ n (t) a stochastic term of 
zero average and autocorrelation function fixed as 



aD 



SnmS(t - t') 



(4) 



and modeling the effect of a thermal bath in contact with 
the system. In the following we will investigate mainly 
deterministic discrete ratchets in absence of the d.c. bias 
7 = and in presence of the a.c. driving of zero average 
([3]). We shall introduce external noise in the system only 
to check that the directed motion is associated to single 
stable attractors, a fact which permits to avoid averaging 
on initial conditions (e.g. initial phase 9 of the driving). 
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FIG. 2: (Color online), (a) Peierls-Nabarro relief for A = 
(continuous line), A = 0.4 (dashed line) and A = 0.9 (dashed- 
dotted line), (b) Peierls-Nabarro energy versus lambda, (c) 
Dependence of the kink width on the asymmetry parameter. 



In Fig. [lja) we show the potential chain in ([!]) for the 
case A = 0.5. Notice that minima of the potential are lo- 

1-V1+8A 2 



4A 



cated at n — S(X), being S(X) — ^- sin 

Since a discrete kink (anti-kink) can be thought as a 
string in space extending from one absolute minimum of 
the potential to an adjacent one, we will assume A S [0, 1] 
to avoid the appearance of relative minima in the poten- 
tial. For large coupling constant K, approximated analyt- 
ical expressions for kinks and anti-kinks can be obtained 
from the corresponding exact solutions of the continuum 
limit derived in |10j (notice that for large k the discrete 
chain is well approximated by the continuous double sine- 
Gordon equation), as 



two discrete kinks obtained numerically for the case of 
symmetric (A = 0) and asymmetric (A = 0.5) potentials. 
Notice that in the symmetric case the kink is symmetric 
around the center [this is clearly seen from the derivative 
profile shown in panel |ljc)] while in the case of asym- 
metric potential the kink is asymmetric and the width is 
a little shrunk. 
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FIG. 3: (Color online), (a) Small oscillations frequency spec- 
trum versus A. The wide band denotes the continuum spec- 
trum, the curve which splits off from the band corresponds to 
the internal mode while the bottom line refers to the trans- 
lational mode (zero mode). The profiles of the internal mode 
and of the translational mode at A = 0.46 are shown in panels 
(b) and (c) respectively. The value of the coupling constant 
k is fixed to 1. 



6(X) 



tan 



sign(A) A B 



A — I — rj sinh 



\ I AB 
2 V l A l 

(5) 

= 2A^2(1 + A), B = 
-/^jpa (the plus and minus 



with A = y/l + 8A 2 , 77 
V /2(4A 2 - 1 + A) and r„ = 
signs refer to the kink and antikink solutions, respec- 
tively). For small values of k the kink solutions must be 
found numerically. In the panel (b) of Fig. [l] we depict 



A. Peierls-Nabarro barrier and depinning 

A main feature of DSR with respect to continuous soli- 
ton ratchets is the presence in the discrete case of the 
Peirls-Nabarro (P-N) barrier which introduces a depin- 
ning threshold for kink transport. To investigate the de- 
pendence of this threshold on the system parameters we 
consider static kink solutions of the unperturbed system 
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solving numerically the nonlinear system of equations 

V'(x n ) + n(2x n - x n+1 - x n -x) = G, n G [-N,N], 

(6) 

with the kink boundary condition xpj — £-iv = +1- 
Starting from a stable kink state we can obtain other 




FIG. 4: (Color online). Small oscillations spectrum versus « 
for A = 0.46. Notice the disappearance of the internal mode 
associated with the oscillations of the kink width for k < 0.20 

(unstable) kink profiles whose center of mass is shifted 
with respect to the stable state, and then find the depen- 
dence of the kink energy on the position of its center 





2000 4000 6000 8000 10000 12000 



x r = 



n=N 

E ■ 

n=-N 



Xn+1 



(7) 



In Fig. |2ja) we show the resulting Peierls-Nabarro (PN) 
relief for three different values of the asymmetry param- 
eter: A = 0, A = 0.4 and A = 0.9. Note that while 
for A = the PN potential possess the inversion center 
symmetry Vpn(X) = Vpn(—X) this is not true for the 
case A ^ 0, due to the asymmetry of the potential in 

0. In all cases the P-N potential is a periodic func- 
tion: Vpn(X + Xq) = Vpn(X), with period X equal to 

1. The difference between the maxima and minima po- 
tential energies, AEpjy, represents the activation energy 
necessary for putting the kink in motion. In panels (b) 
and (c) of Fig. [2] we have also plotted the dependence 
of the activation energy and of the kink width defined 
as [27] 



W K = 1 + 2 



\ 



n=N 

E 

n=—N 



(n - x c y 



(8) 



on the asymmetry parameter A respectively. One can see 
that while the activation energy is an increasing mono- 
tonic function of A, the kink width decreases with A in 
agreement with what observed for the kink profiles in 
Fig. [I] It is worth noting that the above symmetry 
properties of the P-N potential play important role for 



FIG. 5: (Color online), (a) the upper blue line represents 
phase-locked ratchet dynamics of the kink center of mass in 
the absence of noise for A = 0.46 and e = 0.25. The inset 
shows the perfect synchronization of the motion in a smaller 
time scale. When some noise (D = 10" 4 ) is added the di- 
rection of motion is preserved (bottom red line), (b) Chaotic 
ratchet dynamics for A = 0.46 and e = 0.24. Upper blue line 
corresponds to the deterministic case while bottom red line 
corresponds to stochastic dynamics with D = 10 -4 . In all 
cases, a = 0.1, u) = 0.1 and k = 1. 



the existence of directed kink currents in the system. In 
the lowest approximation, indeed, we have that the dy- 
namics of the kink in ^ can be described in terms of an 
effective point-particle equation of the type [B] 



d 2 X{t) 



dt 



, +aX(t) + V' PN (X)+F(t) = 0, 



(9) 



where X(t) and X(t) are collective coordinates denoting 
the position and the velocity of the kink, respectively, 
and F{f) is an effective driving force proportional to F(t) 
in Eq. In order to have directed transport for the 

kink one must break all symmetries that change sign to 
the kink velocity and leave the equation of motion in- 
variant [3 [6] (the presence of such symmetries would im- 
ply the presence in the phase space of kink trajectories 
with opposite velocities, a fact which is obviously incom- 
patible with ratchets). Since the P-N potential in our 
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case does not have the center inversion symmetry and 
since the force is symmetric in time the only symmetry 
which could link trajectories with opposite velocities is 
the one which is indeed an exact symmetry of Eq. |9]) 
in absence of damping. In presence of damping, how- 
ever, the time inversion symmetry is broken and so are 
all the symmetries of Eq. ([9]) which link orbits with op- 
posite velocities. This simple point particle symmetry 
argument shows how important is the damping in our 
system to establish directed transport. It should be re- 
marked, however, that kinks are not exactly point parti- 
cles but extended objects with internal structure which 
also plays an important role for directed transport, as we 
shall discuss in the following. 




FIG. 6: (Color online). Depinning threshold for a dc bias 
7 as a function of the asymmetry parameter A. Parameters: 
a = 0.1, oj = 0.1, £ = 0, k = 1. 



section this prediction is confirmed by simulations cor- 
roborating the relevance of the internal mode in the phe- 
nomenon. 
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B. Internal mode of discrete kinks 

As previously mentioned, discrete kinks of the double 
sine-Gordon system possess an oscillation of the kink pro- 
file (internal mode) even in the pure Hamiltonian limit. 
To investigate the existence of internal modes we have 
linearized the unperturbed system ^ around the static 
kink solution. The numerical spectrum obtained is shown 
in Fig. [3]ja) as a function of A. One can see that, apart 
the zero mode existing for all values of A depicted in 
panel (c), an internal mode emerges for A ^ below the 
phonon band. Its shape, plotted in panel (b), shows a 
spatial asymmetry induced by the asymmetry of the po- 
tential. This mode plays a crucial role in the appearance 
of directed motion through its coupling with the transla- 
tion of the kink. 

It is worth noting that the existence of that internal 
mode also depends on the value of the coupling constant 
K. In Fig. [2] it is shown that for a fixed value of A the 
internal mode disappears for low enough coupling. In 
principle one can predict that no kink net motion must 
exist below that critical coupling. As we will see in next 



FIG. 7: (Color online). Dependence of (v) with respect to A 
for e = 0.1 (a) and e = 0.2 (b). Full circles (crosses) corre- 
spond to phase-locked (chaotic) dynamics. Numbers in brack- 
ets corresponds to the (l,k) pair of Eq. (10 1. In both panels 
a = 0.1, u) = 0.1 and k = 1. 



III. PHASE-LOCKED AND CHAOTIC DSR 

As mentioned before, directed kink motion is expected 
when symmetries of Eq. ([2| that relate kink solutions 
with opposite velocities are broken E]. This is the 
case for the double-sine Gordon potential ([I]) which vi- 
olates the symmetry V(x) — V(—x). In order to ver- 
ify the existence of kink transport and to study the de- 
pendence of the phenomenon on system parameters, we 
have performed numerical integration of Eq. ^ by using 
a fourth-order Runge-Kutta method with constant step 
size At = T/200, being T the period of the external ac 
force. Our system is composed of 400 particles and in 
our simulations static kink was taken as initial condition 
and aperiodic boundary conditions — £-jv = 1 were 
also adopted. 
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Three main dynamical regimes are observed: i) no net 
motion; ii) phase-locked directed motion; iii) chaotic di- 
rected motion. Fig. [5] shows two examples of both kinds 
of movement. In all cases we checked that the phe- 
nomenon does not depends on the initial phase 9 of the 
driving force ([3]), e.g. it is independent on the choice of 
the initial conditions so that one can properly talk about 
ratchet effect. 
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FIG. 8: (Color online). Dynamical regimes of the kink mo- 
tion for different values of the asymmetry parameter A and 
the driving amplitude e. Full circles (crosses) correspond 
to phase-locked (chaotic) dynamics. Blank spaces stand for 
pinned-kink states. Parameters: a = 0.1, ui = 0.1, « = 1. 

In the phase-locked dynamics, the kink travels k sites 
during I periods of the force, and therefore its velocity is 
given by 



In this regime oscillations on the kink profile are perfectly 
synchronized with the motion of the kink center of mass 
so that the kink profile is completely reproduced after 
I periods except for a shift in space. This suggests the 
existence of a coupling between internal and translational 
modes of the kink, similarly as for the continuous case. 

In the chaotic regime the kink jumps backward and 
forward in an unpredictable manner and the directed mo- 
tion can only be observed on a large time scale. In this 
regime the kink motion is of diffusive type and resembles 
the one of a brownian particle. We remark that the direc- 
tionality of the motion in both cases is preserved also in 
presence of a noise term of type Q, as one can see from 
Fig. [5j This is a consequence of the fact that in the phase 
space there is only one attractor corresponding to the net 
motion of the kink (in all investigated cases we have not 
found phenomena of hysteresis signaling the presence of 
more attractors for the same set of parameters). This 
allows to concentrate on the deterministic phenomenon 
and compute the mean kink velocity by averaging on time 
but not on the initial conditions. 

Although in principle the spatial symmetry of our sys- 
tem is broken for any value of A 7^ 0, a distinctive feature 



of discrete systems is the existence of a depinning thresh- 
old above which the transport can occur and below which 
the kink is pinned to a lattice site. We have used to com- 
pute it a dc bias 7 instead of an ac force. The result is 
plotted in Fig. [6] as a function of the asymmetry param- 
eter A. The depinning threshold increases monotonously 
with A in accordance with the rise of the P-N barrier 
shown in Fig. [2] For a fixed amplitude e of the ac force 




FIG. 9: (Color online), (a) Phase-locked (full circles) and 
chaotic (crosses) dynamical regimes of the kink ratchet motion 
in the (e, a) plane. Blank spaces stand for pinned-kink states. 
Other parameters are: A — 0.46, w = 0.1 and k = 1. (b) 
Fraction of phase-locked attractors found in panel (a) for each 
value of the damping. 

the dependence of the mean kink velocity on A is shown in 
Fig. [7] It is interesting to note that even for driving am- 
plitudes above the depinning threshold a minimal value 
of the asymmetry parameter is required for the directed 
motion to start. By further increasing A a rather complex 
scenario appears. Chaotic transport plays a dominant 
role. Notice that the mean drift velocity on chaotic orbits 
can be comparable in some cases with the one achieved 
in the phase-locked regime. Kink dynamics is very sensi- 
tive to parameter values and small changes of them give 
raise to sudden changes of the dynamical regime or can 
even cause a current inversion. This last finding is actu- 
ally a striking result since in the continuous case the sign 
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FIG. 10: (Color online), (a) Phase-locked (full circles) and 
chaotic (crosses) dynamical regimes of the kink ratchet motion 
in the (A, a)-plane. Other parameters are: e = to = 0.1, « = 1. 
(b) Fraction of phase-locked attractors found in panel (a) for 
each value of the damping, (c) Modulus of the averaged mean 
drift kink velocity achieved on the attractors in panel (a) for 
A = 0.46 as a function of the damping. 



of A completely determines the sense of net motion [19] . 
Nevertheless, we have checked that the symmetry prop- 
erty (v(X)) — —{v(—X}) existing in the continuous case 
is preserved in the case of phase-locked dynamics. The 



window of mobility widens of course if we increase e and 
more resonances appear [see Fig. Ujh)] but the scenario 
is qualitatively the same. In Fig.^8] we have shown the 
distribution of regular, chaotic and no transport regimes 
in the (e, A)-plane. We see that whereas the regular kink 
transport occurs mainly at small values of A (e.g. in the 
region A < 0.3, 0.3 < e < 0.5), chaotic directed motion 
becomes dominant at large values of A. Notice also the 
presence of a large pinning region for small values of A 
which extends up to intermediate values of e. 

The influence on the phenomenon of the damping con- 
stant, a, has also been investigated. In Fig. [9]ja) we have 
plotted a diagram with the phase-locked attractors found 
on the (e, a) plane for fixed A = 0.46. We have used in the 
numerical calculations a step of 0.01 in the damping coef- 
ficient and a thinner step of 0.005 in the amplitude of the 
driving force. A band structure emerges as a increases 
and, significantly, no phase-locked attractors are found 
when a — ¥ 0, showing the importance of damping in the 
coupling between the translational an internal modes of 
the kink [HI [20]. In panel (b) of Fig. [9] we show the 
fraction of phase-locked attractors found for each value 
of damping; from it, we can observe a non- monotonous 
behavior of the curve and an higher probability of phase 
locking orbits at intermediate damping. 

In Fig. 10 'a) we report chaotic and phase-locked at- 



tractors found by scanning the (A, a)-plane by keeping 
the rest of parameters fixed. Again, chaotic behavior 
is clearly dominant for low damping while the highest 
probability to find a regular attractor occurs at inter- 
mediate values of the damping, as one can see in panel 
(b) of Fig. " 
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Since in the undamped limit the kink 
mean velocity must go to zero (due to the driving force 
which is symmetric in time) , and since for large damping 
the internal mode mechanism is practically suppressed 
(this strongly reducing the drift velocity), one can ex- 
pect transport not only be more probable but also more 
effective at intermediate values of the damping. This is 



indeed what is shown in panel (c) of Fig. 10 where the 



modulus of the averaged mean drift velocities achieved on 
the regular and chaotic attractors are shown as functions 
of the damping. The fact that at intermediate damp- 
ing the coupling between internal vibration and center of 
mass motion becomes optimal also correlates with similar 
studies performed on continuous [TO] IE] and discrete [55] 
soliton ratchets, as well as, with very recent simulations 
[2 8) of discrete kink ratchets in an asymmetric potential 
of polynomial type free of the Peirls-Nabarro barrier. 

The dependence of the mean kink velocity on the driv- 
ing frequency is characterized by the alternation of chaos 
with a series of phase-locking resonances. This is shown 
in Fig. |ll[ a) where we have depicted the kink mean veloc- 
ity normalized to the driving frequency (v)T = 2t:{v)/uj 
so that for phase-locked dynamics (v)T coincides with 
the ratio of rotation numbers k/l. We do not find the 
monotonic decay observed in the discrete SG model with 
biharmonic forces [55] but, as in that case, kink trans- 
port disappears for high enough frequencies u> > 0.75. 
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tant as k increases. Quite surprisingly, wide regions in 
parameter space in which the kink becomes pinned and 
the net motion suppressed appears for intermediate val- 
ues of the interaction constant. This feature was not 
observed in DSR systems driven by biharmonic forces 
and indicates that the transport may be influenced by 
the local dynamics even if the driving amplitude is above 
the depinning threshold for a fixed value of the coupling 
constant k. This unexpected transport suppression could 
be due to anti-resonances with the external driving or to 
chaos in the kink motion inside the local Peirls-Nabarro 
potential which prevents the escape. The precise mecha- 
nism underlying this phenomenon, however, deserves fur- 
ther investigations. 





FIG. 11: (Color online), (a) Dependence of the normalized 
kink velocity {v)T on the driving frequency for e = 0.1, a = 
0.1, A = 0.46 and « = 1. (b) Phase-locked (full circles) and 
chaotic transport (crosses) in the (A, cj) plane. Blank spaces 
stand for pinned-kink states while the green line represents 
half the frequency of the unperturbed internal mode at the 
same value of A. 



FIG. 12: (Color online). Phase-locked (full circles) and 
chaotic transport (crosses) in the (e, ft) plane. Blank spaces 
stand for pinned-kink states. Parameters are: e = 0.1, a = 
0.1, cj = 0.1, A = 0.46. 



IV. CHAOTIC DSR CONTROL 



Notice that the largest resonant step is achieved around 
uj = 0.5 at the rotation numbers k = — 1, 1 = 1, and that 
at low frequencies (e.g. uj < 0.2) the largest drift veloc- 
ity is achieved in correspondence of chaotic trajectories. 
Also notice in panel (b) of Fig. 11 that around to = 0.5 
the phase-locked attractors cover a large region of the 
(A — cj)-plane which partially intersects the first subhar- 
monic curve of the unperturbed internal mode frequency, 
this indicating a coupling between the driving frequency 
and the internal mode Qi « 2co. 

Lastly, we have investigated the dependence of the 
ratchet phenomenon on the interaction constant k. In 
our model this parameter measures the discreteness of 
the system, corresponding k — > oo to the continuum 
limit. As expected (see Fig. 12) high discreteness pre- 



vents propagation of the kink, and the depinning ampli- 
tude of the driving decreases as n increases. Just above 
the depinning amplitude chaotic transport predominates 
while regular phase-locked motion becomes more impor- 



In contrast to continuous soliton ratchet systems 
which, except for the case of strongly localized inho- 
mogeneities [35], are very robust against chaos, discrete 
soliton ratchets exhibit chaos particularly in the case of 
asymmetric potential we are dealing with. We have in- 
vestigated conditions which allow to reduce the chaotic 
orbits in favor of the phase locked dynamics. To this re- 
gard we have added to the sinusoidal driving ([3| a weak 
subharmonic signal 77 cos(aji/2) in phase with it. 

Stabilization effects induced by external weak subhar- 
monic signals are known to be a general mechanism for 
suppressing chaos 30 , both in point-particle systems 
[31] and in infinite degree of freedom systems described 
by partial differential equations [35]. The small subhar- 
monic driving is used not to create new ratchet-like or- 
bits but to stabilize regular orbits over larger ranges of 
the system parameters. 

In our study we have chosen e as bifurcation param- 
eter since, in principle, in a possible experimental im- 



9 




2B 0.23 0.232 0.234 0.236 0.238 0.24 



0.14 
0.12 
0. 

i- 0.08 

T3 

O 

E 

=- 0.06 
0.04 
0.02 




(a) 


% 






*5t. ■ ' ' 








5 — \ 
















; I ' 






% ■ 





0.2 


0.4 


0.6 


o.a 



0.12 
0.1 
*- 0.0B 

O 

E 

0.06 




i)'"2 



28 0.23 0.232 0.234 0.236 0.238 0.24 




FIG. 13: (Color online). Bifurcation diagrams for y mod 1 
with respect to e for 77 = (a) and r\ — 10 -4 (b). Variable y 
mod 1 is obtained taken the position of the kink center each 
period of the driving after a transient. Other parameters are: 
A = 0.4, a = 0.1 and cj = 0.2. 



FIG. 14: (Color online), (a) bifurcation diagram for y mod 1 
with respect to parameter m of the elliptic function (111 and 



e = 0.232. The amplitude e' has been adjusted according to 
(13 1, (b) bifurcation diagram as a function of e for m — 0.7. 
In both panels parameters are the same as in Fig. |13| 



plcmcntation of the system, e would be easier to control 
and change than A for instance. For simplicity we have 
focused our study on an interval of e where a destabi- 
lization of phase locking orbits occurs. In Fig. [i"3^a) we 
have plotted, for r\ = 0, a variable y mod 1 versus e, where 
y is defined as the stroboscopic map of X c with period 
T = 2tt/uj (i.e. X c taken each period T after a transient). 
One can observe phase-locked orbits represented by con- 
stant values of y mod 1 for a given e. Transition to chaos 
occurs via a cascade of period doubling bifurcations. In 
this graph, for the sake of continuity of the variable y, we 
have used as initial condition for each step the solution 
obtained at the previous one. Another stroboscopic map 
for the same interval and 77 = 1CP 4 is shown in panel (b) 
of Fig. [13] We see that in this case the subharmonic 
signal, although significantly reduces the dispersion in y 
in the chaotic region, it does not give the full suppression 
of the chaotic region. This suggests that other harmonic 
signals may also be involved in the full chaos suppression. 

To find the optimal shape of the small signal for chaos 
suppression, it is convenient to use a more general driving 
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FIG. 15: (Color online). Bifurcation diagram as a function of 
e obtained for m = 0.7 when the elliptic function is replaced 
with the first two terms of its Fourier expansion. Parameters 
are the same as in Fig. [13] 
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in the form of a Jacobi elliptic function 
F(t) =£cnl ojt, to 



(11) 



with e' the amplitude of the driving, < to < 1 is the 
elliptic modulus and K is the complete first order ellip- 
tic integral K{m) = f* /2 (1 - msin 2 6)- 1 / 2 d9 33 . This 
driving allows us to change the periodicity of the sig- 
nal continuously by means of the modulus m used as a 
free parameter to optimize the shape. An elliptic driv- 
ing, apart from being periodic, satisfying the symmetry 
F(t) = —F(t + T/2) and reducing to a sinusoidal driv- 
ing ([3]) in the limit m — 0, has the advantage of being 
expressed in terms of sinusoidal functions as 



F(t) = 



2tt 



00 



n+1/2 



m i/2 K 1 

n=0 



- cos[(2n + l)wt] , (12) 



with the coefficient of the series depending on the param- 
eter q = e -T*K(i-m)/K{m) an( j orn y qq}^ multiples 

of the fundamental frequency entering the sum. Similar 
drivings were also used to suppress chaotic dynamics of 
Bose-Einstein condensates loaded into a moving optical 
lattice [34] , 

The above driving was used to stabilize, for exam- 
ple, the chaotic region of our original system around 
e = 0.232 in Fig. 13 ^a). By replacing the sinusoidal sub- 



harmonic driving with the elliptic one (11) and starting 



from to = 0, we have increased m looking for a regular- 
ization route to a phase-locking dynamics. In order to 
properly compare with the sinusoidal case (m = 0) we 
have adjusted, while increasing to, the amplitude e' so as 
to keep the first Fourier coefficient of the expansion ( 12 ) 



to be e, i.e. we performed the regularization by following 
the path 



e'(m) 



if 771 — 

if0<m<l 



(13) 



27iy/2 



The results are depicted in Fig. 14 ^a). We see that as 
to is increased a progressive regularization occurs. The 
system undergoes an inverse period doubling sequence 
of bifurcations until it reaches a periodic orbit around 
777 = 0.7. We have chosen this value to compute now a 
new stroboscopic map with e as bifurcation parameter 
which it is shown in panel (b) of Fig. 14 Comparison 



with Fig. 13 '&) reveals the complete regularization of the 
system for e > 0.2316. 

An elliptic driving in practical applications could be 
difficult to implement. We observe, however, that the 



main contributions in the series expansion (12 1 come 



from the first and third harmonics. Thus, for instance, 
for m = 0.7 the coefficients of the third and fifth har- 
monics are 0.0802e' and 0.00603e' respectively, being e' 
the amplitude of the first harmonic. In Fig. [l5]we have 
depicted the resulting stroboscopic map when the ellip- 
tic function (11) is substituted for the first two terms of 



its fourier expansion. Note that the addition of a third 
harmonic to the fundamental frequency is indeed enough 
to suppress chaos. Therefore a reasonable option in a 
possible experimental set-up to control chaos would be 
to add to the ac driving a third harmonic conveniently 
adjusted. 



V. CONCLUSIONS 

We have studied the ratchet phenomena induced by 
asymmetric potentials on discrete topological solitons 
(kinks and antikinks) of the discrete damped double sine- 
Gordon equation in the presence of an harmonic driving. 
The dependence of the mean velocity of the kink on the 
asymmetry parameter, damping, amplitude of the driv- 
ing, frequency of the driving and on the coupling con- 
stant, has been extensively investigated. We have shown 
that in contrast with DSR driven by asymmetric forces of 
zero mean |22j . kink transport can be quite effective not 
only on regular (phase-locking) orbits, as it is usual for 
soliton ratchets [10], but also on chaotic trajectories. In 
particular, we have showed that although the mean kink 
velocity in the chaotic region is typically smaller than the 
one achieved on phase-locked orbits, there are parameter 
values for which the chaotic transport becomes more effi- 
cient than the one in the phase locking regime (this case 
is observed, for example, in Fig. 11 'a) for small values of 
the driving frequencies, e.g. for ui < 0.2). 

Our investigation of the dependence of the kink trans- 
port on the damping constant shows that a phase lock- 
ing regime and therefore an efficient transport is more 
probable at intermediate values of damping. This fact 
correlates with previous investigations and explanations 
in terms of the internal mode mechanism. 

The mechanism of the kink transport also depends 
strongly on the coupling constant. We find that for large 
values of the coupling constant the kink transport is dom- 
inated by the regular phase locked motion while the diffu- 
sive chaotic transport becomes important for small values 
of k. Quite surprisingly, we have also found that for in- 
termediate values of k a wide region of pinned kink states 
can exist. 

The ratchet dynamics in the chaotic regime has been 
investigated in terms of a reduced map for the center of 
the kink which allows us to show that the onset of chaos 
occurs via period doubling bifurcations. The possibility 
to control chaotic ratchet dynamics has also been inves- 
tigated by means of an additional sub-harmonic driving 
which does not break the time symmetry of the system. 
In particular we have shown that the chaotic region can 
be strongly reduced in favor of periodic phase-locked or- 
bits even for very small strengths of the subharmonic 
signal in comparison with the fundamental driving. To 
fully suppress chaos we have used a more general driving 
in the form of an elliptic function with the modulus taken 
as a free parameter to modulate the signal and to find 
the optimal regularization route for the system dynam- 
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ics. For the region of parameters investigated, we have 
found that in any practical context this would lead to 
phase-locking ratchet transport and to chaos suppression 
simply by adding a third harmonic conveniently adjusted 
to the fundamental driving. 

Finally, we remark that the features of DSR investi- 
gated in this paper may be of general validity also for 
other systems undergoing structural changes which in- 
duce an asymmetry in the potential energy. 
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